Quantum Particle Production 
at Sudden Singularities 



oo : 
O . 

O ■ John D. Barrow^, Antonio B. Batista^, 

^ ! Julio C Fabris^'^ and Stephane Houndjo^ 

^ ' ^DAMTP, Centre for Mathematical Sciences, 

\-^ • Cambridge University, Wilberforce Road, 

^ '. Cambridge CBS OWA, UK 

^ ! ^Departamento de Fisica, 

I Universidade Federal do Espirito Santo, 

^' CEP 29060-900 Vitoria, Espirito Santo, Brasil 

O^. ^ GreCO, Institut d'Astrophysique de Paris - lAP, 

;-i ! 98bis, Boulevard Arago, 75014 Paris, France 

Cn ■ Abstract 

> ■ 

^f-\ ' We investigate the effects of quantum particle production on a classi- 

1/^ ^ cal sudden singularity occurring at finite time in a Friedmann universe. 

^SJ ■ We use an exact solution to describe an initially radiation-dominated uni- 

^^ ' verse that evolves into a sudden singularity at finite time. We calculate 

r ^ , the density of created particles exactly and find that it is generally much 

f^ ' smaller than the classical background density and pressure which produce 

QQ ' the sudden singularity. We conclude that, in the example studied, the 

f^ , quantum particle production does not lead to the avoidance or modifica- 

ll" ' tion to the sudden future singularity. We argue that the effects of small 

. 5^ I residual anisotropics in the expansion will not change these results and 

^^ . show how they can be related to studies of classical particle production 

$_( ' using a bulk viscosity. We conclude that we do not expect to see signif- 

5t ! leant observable effects from local sudden singularities on our past light 
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1 Introduction 



In ref. [1] it was shown that classical general relativistic Friedmann universes 
allow finite-time singularities to occur in which the scale factor, a(t), its time 
derivative, a, and the density, p, remain finite whilst a singularity occurs in the 



fluid pressure, p -^ +00, and the expansion acceleration, a -^ 00. Remarkably, 
the strong energy condition continues to hold, p + 3p > 0. Analogous solutions 
are possible in which the singularity can occur only in arbitrarily high derivatives 
of a(t),[2]. This behaviour occurs independently of the 3-curvature of the uni- 
verse and can prevent closed Friedmann universes that obey the strong energy 
condition from recoUapsing [3]. Subsequently, a number of studies have been 
carried out which generalise these results to different cosmologies and theories 
of gravity [4] , classify the other types of future singularity that can arise during 
the expansion of the universe [5] , and explore some observational constraints on 
their possible future occurrence in our Universe [6] . 

In this paper we extend these studies to consider some quantum implica- 
tions of a sudden singularity. Specifically, we want to know if quantum particle 
production can dominate over the classical background density on approach to 
a sudden singularity and stop it occurring or modify its properties, as can be 
the case for the Big Rip future singularities [7, 8] . The results of such a study 
are also of interest for observational probes of finite-time singularities. In an 
inhomogeneous universe it would be possible for us in principle to observe a 
sudden singularity on our past light cone. What might we see? If quantum ef- 
fects produce profuse particle production then there might be observable effects 
from local sudden singularities. In this respect, sudden singularities are also of 
interest with regard to proposed measures of gravitational entropy associated 
with invariants of the Weyl curvature [9] . The Weyl invariant will not diverge on 
approach to a sudden singularity (and there is no geodesic incompleteness [10]) 
and so it may represent a part of a soft future boundary of the universe with 
low gravitational entropy - which could be as close as 8.7 Myr to the future [6]. 
The effects of loop quantum gravity have been studied in cosmologies exhibiting 
classical sudden singularities and they can remove the sudden singularity un- 
der certain particular conditions [11], and there is a close relationship between 
sudden singularities and the behaviour of Friedmann universes containing bulk 
viscous stresses [12]. 

In order to provide some insights into these issues we will construct a simple 
exact classical Friedmann cosmological model with a future sudden singular- 
ity that follows an era of radiation domination in which there is no quantum 
production of masslcss particles. We will calculate the quantum production of 
massless particles on approach to the sudden singularity where there is a de- 
parture from conformal invariance. For simplicity, wc will ignore any period 
of late-time acceleration in the universe although this more realistic detail can 
easily be incorporated. We will then discuss various elaborations of this model 
and show why we do not expect small deviations from isotropy and homogeneity 
to alter our conclusions. 



2 A Suddenly Singular Cosmology 

We will employ the simplest example of a spatially-flat isotropic Friedmann 
universe with a sudden singularity, presented in [1] . The cosmological equations 



for scale factor are: 
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We take the following solution for the scale factor and its first and second 
derivatives: 
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Imposing < g < 1 and 1 < n < 2, we sec that a -^ a^Q, a -^ iy°-^o- ) ^^^ 
a — > oo, as i -^ tso- Hence, the scale factor and its first derivative remain finite, 
implying that the density remains also finite, while the second derivative of the 
scale factor, and consequently the pressure, diverges. We assume the condition 
< g < 1 in order to have a decelerating universe as i — ^ 0, but this condition 
could be relaxed without affecting our results. 

There is no exact solution for the Klein-Gordon equation for the form (4) 
for the scale factor. Therefore, we consider some simplifications. We divide the 
entire evolution of the universe into two phases: one which characterizes the 
primordial phase, t ^ 0, and other which characterizes the "singular" phase, 
i —^ tso- For the primordial phase, we will use the radiation-dominated phase 
of the standard model (i.e. q — 1/2), for which it is possible to solve the Klein- 
Gordon equation so that the solution naturally contains the structure of the 
vacuum state of quantum fields. 

Let us consider the asymptotic behaviours of the solution (4): 



• Primordial radiation phase (i ^ 0): 



«sO - 1^1/2^ (7) 



1/2 



(a,o - l)_^_i/2^ (8) 



2t'/' 



Singular phase {t -^ tso)' 
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We introduce the conformal time, r], defined by adi] = dt and re-express 
the scale factor during the radiative phase as a = rt^, where the constant 
r = "°°~^ , so ?7 = 2^2. This leads finally to the following expressions: 

• Primordial radiation phase (77 — > 0): 

a == aoT], (11) 

a' = ao- (12) 



Singular phase (ry -^ Vso)- 



a = aso, (13) 
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In the above expressions, ao = ^ and ' = d/drj. The transition to the 
singular phase now occurs at the conformal time rj = tjq and wc have 

apiVo) = asirjo), a'piij^) ^ a'^ivo), (15) 

where the subscripts p and s denote primordial and singular phases, respectively. 
This implies that 
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The matching conditions imply 

ao = -^, aso = 2^/t7o, (17) 

^0 

and the solutions for the two phases arc: 

1. Primordial radiation phase: 

a == flso — , a = ao = ; (18) 

Vo Vo 

2. Singular phase: 

a = aso, a = ao = . (19) 

^0 

In order to construct the simplified model of approach to a sudden singular- 
ity, we will consider 770, ry^ and aso to be independent free parameters. 



3 The Klein-Gordon equation 

In a spatially-flat Friedinann universe, the Klein-Gordon equation for a massless 
field, 6, is 
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where we have re-inserted the light velocity, c, and if k is the wave number of 
the Fourier decomposition: 
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For simplicity, we have omitted the Fourier index in the function (p. Now, we 
scale by H^, where Hq is the Hubble parameter at the moment of the transition, 
rjQ. Moreover, we define k = kin, where Ih is the Hubble radius at tjq. Note 
also that the scalar field has dimensions of (length) ^/^. 

Now we consider the Klein-Gordon equation for the two phases defined 
above. For simplicity, we will omit the bars, setting fc — > fc, and Horj -^ rj 
(a dimensionless time parameter). Note that with this parametrization, t^q = 1 
(since i^o^o = 1 in the old variables). 

3.1 Primordial phase 

In this case, 0(77) = a^rj. Hence, 
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The solution of this equation is 
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where Hi ' (x) are the Hankel functions of the first and second kind. 
The Hankel functions are defined as follows: 



H[]lix) = J,/2ix)+iN,/2ixy,H[f^ix) - Ji/2(x) -j7Vi/2(x), (24) 

where Jy{x) and N^{x) are the usual Bessel and Neumann functions. Using the 
fact that 
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n x 
the solution for the scalar field can be re-written as 
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We choose 



Ci-^^A/^^p, C2=0, (27) 



where the factor \/ ^Ipi has been chosen due to the dimension of the scalar field 
for later convenience, and obtain the typical behaviour of a normalized quantum 
vacuum state with a time factor added: 

MV,^)-^^^^^. (28) 

3.2 Singular phase 

In the singular phase, the Klein-Gordon equation reduces to 

</)" + 2(1)' + k^(j) = 0. (29) 

The solutions are 

^^A+eP+^ + A^eP-", (30) 

where 

p± = -1± V^^-l- (31) 

The final solution can be written as 



^kiv^x) = e "^ 



j\ gi{i^V-k-x) _|_ ^^g-i(wij+fe-K) 



(32) 



where u> = \/k'^ — 1, and so the field propagates only if /c > Hq. We have also 
incorporated the factor \j\lpi in the definition of the constants A±, to make 
them dimensionless. 

3.3 The potential 

One way to visualizing the overall solution is to redefine the scalar field so that 

<^=-, (33) 

a 

and then Klein-Gordon equation takes the form 

A" + [fc2 - V{a)\\ = 0, (34) 

where Via) — ^- Hence, there is propagation whenever fc^ > V(a), and the 
field decays when fc^ < V{a). That is, we have a quantum mechanical problem 
of particle propagation in a potential barrier. 

For the radiative universe, a (X r], leading to V{a) = in the primordial 
phase: all modes propagate [13]. In the singular phase, however, we have a 



more complicated situation, since we must evaluate the second derivative of the 
scale factor near the singularity. We find, 
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Since 



we find that the potential for this phase is 



n«) = - = -n(n-l)^^, ^Bl (37) 

This confirms that there is no propagation for fc < Hq. 

4 Matching the solutions 

Using the solutions of the Klein-Gordon equation in the two phases, we match 
the fields and their first derivative at 77 = ryg. This gives the following expres- 
sions: 

A+e*'^ + A_e-^'^|, (38) 



{-\ + ik)\\^e'^ = A+(-l + iuj)e^"^-^'> -A^(l + iuj)e-^"^+^\ (39) 
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Taking the combination {{1 + iuj) x (38) + (39), and later ((l-w)x(38) + (39), 
we find: 

5 Energy density of created particles 

The energy density of the created particles is given by (see [14]), 

trrS foo r ^ 

Ps = j^^J^ dfcfc2|<^'^0'; + fc20^<^*j. (42) 

This expression is obtained by computing the vacuum expectation value of the 
Hamiltonian for a quantized massless scalar field. It is completely equivalent to 
the alternative derivation of the energy density using the Bogoliobov coefficients 
[7, 15]. The pressure associated to the created particles is given by [14]: 

'^^o / dkk^^',r,-^^,rA- (43) 
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Now, using (32) and (40,41), wc find the following expression for the energy: 
3hH^ Liey f^^ „ 4fc 
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cos[v/A:2 - ly] + ^k^ - 1 sin[v/A:2 - ly] I (44) 



where j/ = 2(1 — ?]). The transition occurs at rj — rjQ — I {y — 0), noting the 
re-scaling. The upper limit of integration was set by fc^, an ultraviolet cut-off 
that must be identified with the Planck wavenumber. The reason for this is 
that we expect the Klein-Gordon equation may not retain its simple form (20) 
at energies higher than the Planck scale, where quantum gravity enters. There 
are many studies of the modification of the usual dispersion relation for a scalar 
field in the context of black hole thermodynamics and cosmological perturba- 
tions of quantum origin, two situations plagued by transplanckian frequencies 
- see [16] and references therein. The modification of the dispersion relation is 
usually treated by introducing a decreasing exponential term, which leads to 
a very important suppression in the integration in the transplanckian region. 
The adoption of this procedure is equivalent to stopping the integration near 
the Planck frequency. Since the extrapolation to the transplanckian regime is 
very speculative, we will ignore this transplanckian problem and adopt a more 
conservative regularisation procedure [15]. We will return to this problem later. 
Let us choose the scaling so that Uso = 1. The above expressions can be 
then be rewritten as 

- cos[\/fc2 - ly] + Vk^ - lsm[^/k^ - ly] I , (45) 

where Pq is the background cosmological density at the time of the transition. 
This integral can be solved exactly and we have 

p, - 2^-^^ey\kt + kl + Hkl-l]-2Ci[-^i^~ly] 
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In this expression, Ci[x\ is the cosine integral function and Chi[x\ is the hyper- 
bolic cosine integral function. 

The general behaviour indicates that the energy density of created particles 
decreases, see the bottom left graphic in figure 1. This expression is rigorously 
valid only after the transition. The initial condition is the number of particles 
created during the first (radiative) phase. The plot is made in units of this initial 
number, which is, at best, small. From this we can conclude that the sudden 
singularity is robust against particle production due to quantum effects, since 



the energy density of created particles is generally much smaller than the energy 
density of the background, and goes quickly to zero as the future singularity is 
approached. This result also confirms the self-consistency of our calculation 
using the Klein-Gordon equation to describe the evolution of a quantum field 
in the metric of a classical background cosmology [15]. 

If we now take the limit km —^ oo, keeping the Klein-Gordon equation in its 
traditional form (ignoring the transplanckian problems) , we find that there is a 
quartic, a quadratic and a logarithmically divergent term. Moreover, the last 
term in the integral (44) may only ultimately have a meaning if it is treated as 
a distribution. Such divergencies may be removed using regularization of the 
energy-momentum tensor. The general result depends strongly on the back- 
ground, which in our case is quite non-trivial. Even if we consider the simplified 
scenario with a radiative phase preceding a singular phase, the regularization 
could be implemented in the radiation-dominated phase (which is a known result 
[15]), but the existence of a discontinuity in the second derivative, due to the 
matching conditions we employ, makes the application of the usual expression 
neither direct nor simple. But, we can proceed in a more straightforward way. 
In the singular phase, we we are evaluating the creation of particles and the scale 
factor becomes constant. This makes it secure to use a covariant subtraction of 
the infinities [17], for the density and for the pressure, 

where < T^i, >reg stands for the regularized energy-momentum tensor, < 
T^u >0i for the energy-momentum tensor evaluated using the expressions above, 
and < T^u >div is its corresponding divergent part. The divergent parts are 
represented by the first term in (44). The final expression is given by (see the 
Appendix) 
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It is very important to stress that the Minkowski limit, obtained by imposing 
Pq —f 0, leads to a null result, the same as we would obtain if we had computed 
the vacuum expectation value in Minkowski space-time and subtracted the di- 
vergencies. At same time, the resulting energy-momentum is conserved. This 
confirms the consistency of the procedure employed above. Note that, in the 
present case, the equation of state does not coincide with the classical equa- 
tion of state of a massless scalar field in a FRW background, which is that of 
a stiff matter fluid {p = p). In general, the quantization of a classical system 
may change the classical equation of state. The general form of the regularized 
energy density of created particles is exhibited in figure 1 , showing a decreasing 
behavior as before. 

It is possible to obtain at least some information about the effects of using a 
simplified model, where the evolution of the universe is described by two phases, 
by integrating numerically the exact Klein-Gordon equation using (4,5). The 
results are of course plagued by the problem of the ultraviolet divergence. But, 
we can stop the integration at a high frequency and compare the result with the 



expression (46), or even with (48). The initial condition is given by the same 
vacuum state that we used above. The results for the background and for the 
energy density of particles created are compared in figure 1 , where we have also 
inserted the expressions for the energy of the produced particles using the Planck 
frequency cut-off (46) and the regularized energy-momentum tensor (48) as well 
as the numerical computation. Of course, the regularized expression (48) fits the 
numerical integration only in a general form; this is natural since a divergent 
contribution has been extracted to obtain (48). Yet they agree in the sense 
that both predict a decreasing energy of produced particles as the singularity 
is approached. This is a quite general behavior that remains even if the free 
parameters, like Ogo and Hq, are varied. As is also expected, the numerical 
fitting to the simplified analytical expression is better when the duration of 
the singular phase is small compared to that of the preceding radiative phase. 
Moreover, the numerical integration reveals the effect of the sharp transition 
(which is not displayed in the figures): in (46,48) there is a divergence in the 
energy at y = 0, which does not appear in the numerical computation. This is 
an effect of the discontinuity in the second derivative, which is clearer when the 
Klein-Gordon equation is written as in (34). 

To obtain more details of the result described above, we can use the expres- 
sion for the particle production exhibited in reference [15]. The energy density 

can be written as 

fir r°° 

Nkkd^k, (49) 



where Nk is the particle occupation number. Comparing with (45), and retain- 
ing only the relevant terms after rcgularization, the particle occupation is then 
given by 



Nk oc ^^ , 1) 1 (~ cos[Vk^ - ly] + Vk^-l sin[v/fc2 -ly][. (50) 

From this expression we can verify that the particle production decreases as 
the sudden singularity is approached, in contrast to what happens in the case 
of approach to a 'big rip' future singularity [18] (where a — > oo,p ^ oo and 
IpI ^- oo as i ^ ts). Moreover, light particles are produced more copiously than 
heavy particles, as would be expected. These behaviors are shown in Figure 2. 



We can try to compare our results with those obtained in reference [8] . How- 
ever, the context is quite different to that of this paper since the authors of [8] 
have considered an ensemble of conformal quantum fields, (whereas we consider 
only a non-conformal, massless scalar field), generating a trace anomaly by us- 
ing the effective-action approach where gravity is modified by requiring that the 
quantum fields must be renormalisable on a given metric background. They find 
that the sudden singularity can be modified by quantum effects. The energy 
conditions could be violated and a big-rip singularity ensue. In our calculation, 
these quantum effects are inoperative. A possible reason for these differences is 
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the expression used in ref. [8] for the conformal anomaly: it was derived near an 
initial strong curvature singularity (a -^ 0, d/a ^-oOjp^ooasi— >0). Since 
the scale factor becomes constant near the sudden singularity (a,d,p all finite 
but a ^ oo as i — > is), it is possible that the expression used for the conformal 
anomaly does not apply unchanged from its form in a scenario with a dynamic 
scale factor near an initial singularity where the density diverges. The quantum 
effects considered in reference [8] have their natural domain of applicability in 
the early universe, near the big bang rather than at a late-time sudden singu- 
larity where geodesies are undisturbed. At a future sudden singularity there 
is a curvature singularity which can in principle justify the employment of a 
conformal anomaly for the sudden singularity. However, the sudden singularity 
has many features which distinguish it from the initial big bang singularity: the 
divergence appears only in the second derivative of the scale factor (or, equiv- 
alently, in the pressure) and the expansion rate and the density remain finite. 
Hence, the validity of such an extrapolation to the sudden singularity remains to 
be proven. In particular, if particle production effects from other type of fields, 
like massive or conformal scalar fields, or even vectorial and spinorial quantum 
fields, produce such a large back-reaction that a sudden singularity is changed 
into a big rip singularity then it is difficult to describe that self-consistently in 
terms of effective quantum stresses on a fixed background because the back- 
ground expansion is strongly perturbed. 



5.1 Classical analogues 

The quantum particle production effects in an isotropic universe can be viewed 
as arising because of the presence of an effective bulk viscosity, ^, in the energy- 
momentum tensor [19]. The presence of a classical bulk viscosity leaves the 
Friedmann eq. (1) unchanged and modifies eq. (3) to 

p + 3H{p + p) = 9H^^, (51) 

where ^{p) — ap"^ on the right-hand side represents a classical analogue of the 
non-equilibrium particle creation effects produced by the viscous stress [12]. On 
approach to a sudden singularity, a -^ a^o, H -^ HgQ and p -^ p^Q, while p and 
\p\ —* oo, so we see that ^ —^ ^^g = ^{Pso) and the classical particle production 
term 9H^^gQ both approach constants as i ^- tso, so to leading order 

p + 3Hp -^ 9H^^ -^ constant, 

and the production effects are classically negligible. 

It is also interesting to ask if the presence of anisotropics in the cosmological 
expansion rate could significantly change our results. In the study of quantum 
effects on approach to a curvature singularity in the early universe we are fa- 
miliar with the strong effects of quantum particle production, which can bring 
about significant isotropisation of the expansion, both directly and as a result of 
the anisotropic red and blue-shifting of the created particles after they become 
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collisionless [20] . These effects are driven by the fast divergence of the dominant 
shear anisotropy energy density (ct^ oc a~^) at small a. We could ask whether 
any possible amplification of small anisotropies on approach to a sudden singu- 
larity could lead to major quantum effects. However, we can easily see that such 
dominant effects will not occur on approach to a sudden singularity as a ^ as 
when t ^ ts- In this situation the anisotropy energy density also approaches 
a constant value, p{ts) oc aj^ which will be far smaller than the ambient den- 
sity of dust or radiation, and the associated density of created particles will be 
smaller still. There will just be small changes to the asymptotic energy density 
and expansion rate in the limit a ^ as, as eq. (1) is modified to 

3Hs = SttG/Sso + o-sO' 

with all three terms equal to constants. The dissipative effects of the created 
particles are analogous to the presence of a shear viscous stress, proportional to 
the shear and its effects also remain small as i ^ iso- 

6 Discussion 

We have considered the quantum particle production that would occur in an 
isotropic and homogeneous universe on approach to sudden singularity, where 
a, a and p remain finite while a and p diverge. If significant, such quantum effects 
could modify or remove a sudden singularity at finite time. We have set up a 
simple exactly soluble model in which an early radiation-dominated universe 
evolves towards a sudden singularity at finite time. We compute the quantum 
particle production as the singularity is approached and show that its effects 
remain negligible with respect to the classical background pressure and density 
all the way into the singularity. We compare our results to other discussions 
of quantum effects in the literature; we argue that any effects created by the 
presence of small anisotropies will be negligible and show how a simple classical 
description of the quantum particle production as an effective bulk viscosity in 
a Friedmann universe gives a similar outcome. 
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Appendix 

The energy density and the pressure are given by (neglecting muhiphcatively 
unimportant terms), 

p, = — e^ / -^\{2k-l)k^ -cos{ujy)+uii\n{uy)\dk , (52) 

Ps - ^eM -^U2k-l) — -il--k^jcos{ujy)+u;sm{Loy)\dk (53) 

In these expressions, w = y/k'^ — 1 and y — 2(1 — rj) where rj is the conformal 
time (which is proportional to the cosmic time t when a ~ constant). Now, let 
us consider the conservation law, 

p' + .i-{p + p)^0 . (54) 

The primes mean derivative with respect to rj. This may be rewritten, after 
redefining the time variable using the (constant) Hubble parameter, in terms of 
the derivative with respect to y: 

-2p + 3(p + p) = . (55) 

It is easy to see that the expressions (46,52) satisfy (55). It is important to note 
that, even if a is constant, we must derive it since d (equivalently, a') is not 
zero, being also equal to a constant. More importantly, if we rewrite (46,52) as 

Ps == Psl + Ps2 , (56) 

Ps = Psl + Ps2 , (57) 



where 



1 r°° I. 

p,i = ^eM —{2k-l)k^dk 



(58) 



1 f°° k 
Ps2 = ^e^ / -^{ -cos{ujy)+ujsm{u;y)}dk , (59) 



a^ Jq uj^ 



Psl 



1 /-oo I, 

-ey / -^{2k~~l)k^dk , (60) 

Ps2 = — e^ / —<-il--k'^jcos{ujy)+ujs[n{ujy)[dk , (61) 

then the pairs {Psi,Ps2) and {pgi,ps2) satisfy (55) separately. 

Now we can easily show that Pg2 is finite. In fact, it can be written as 



Id f°° kdk 



Ps2 



a^ dy Jq ijp' 
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e^ cos(a;j/) 



(62) 



The integral can now be performed, leading to 



1 d 
a^ dy 



Ps2 = -^-(e'CM-y)] , (63) 



where we have exploited the fact that the cosine function is even in its argument. 
This function is finite, except at y = 1, an effect of the sharp transition of the 
second derivative in the simplified model which does not appear in the numerical 
integration, where the background is smooth. 

Now we can obtain the pressure Ps2 ■ The result is 



1 e^ f „, , , cosht/ cosht/ sinhw 

CM{-y) + 3 ^ - 2^^ + 2 '^ 

5a^ I y y^ y 



Ps2^T;-\CM{-y) + 3 -1-2^^ + 2^-^} , (64) 



which goes to zero asymptotically. 

Hence, after subtracting the divergent parts Pg^ and psi, we obtain finite ex- 
pressions for the energy and for the pressure, which have the ordinary Minkowski 
space-time limit satisfying the covariant energy-momentum conservation law. 

The regularization procedure is employed here is inspired in the n-wave 
technique [19] (see also [23]). Here, its application is very simple due to the 
natural splitting of the energy density and pressure into a divergent and a finite 
parts, which separately satisfy the conservation law. 
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Figure 1 : These figures show the evolution of the background scale factor (upper left) of the 
simplified (dashed line) and of the exact (continuous line) models, the behavior of the energy 
of the produced particles using the numerical results (upper right), the truncated (bottom 
left) and the regularized (bottom right) expressions. We used arbitrary normalization. In the 
figures the energies are obtained with different time coordinates (proper time for the numerical 
integration and conformal time for the analytical models) but the future direction is always 
in the increasing value of the coordinate. 
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Figure 2: The density of particles created with a given momentum k as function of time 
(left) and the density of particles created as a function of the momentum k for a fixed time 77 
(right), where kg is a reference momentum (a value feo = 2 is chosen arbitrarily) and Nj^o is 
the corresponding particle density. 
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